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We examine the effect of large extra dimensions on vacuum decay in the Randall- 
Sundrum (RS) braneworld paradigm. We assume the scalar field is confined to the brane, 
and compute the probability for forming an “anti-de Sitter” (AdS) bubble inside a crit- 
ical flat RS brane. We present the first full numerical solutions for the brane instanton 
considering two test potentials for the scalar field. We explore the geometrical impact 
of thin and thick bubble walls, and compute the instanton action in a range of cases. 
We conclude by commenting on a more physically realistic potential relevant for the 
Standard Model Higgs. For bubbles with large backreaction, the extra dimension has 
a dramatic effect on the tunnelling rate, however, for the weakly backreacting bubbles 
more relevant for realistic Standard Model potentials, the extra dimension has little 
impact. 


Keywords: Vacuum decay; bubble nucleation; gravitational instantons. 


1. Introduction 


When Coleman and de Lucci! pioneered the study of vacuum decay in curved 
spacetime, they described the possibility as “the ultimate ecological catastrophe.” 
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Whilst the comment was somewhat tongue in cheek, the recent measurement of 
the Higgs masg2I and the realisation that the Standard Model Higgs field could 
well be in a metastable state#IZ2] has brought the catastrophe a little closer to 
reality! Fortunately, the timescale for decay according to Coleman et al DEEA (see 
also Ref. [I5), is sufficiently large that we would seem not be troubled, except that 
the Coleman results are computed in a highly symmetric background. Recent work 
by two of ud!8l20 has argued that taking into account inhomogeneities such as 
primordial black holes can dramatically shorten the lifetime of the false vacuum 
(see also Refs. BIļand [22]for early work, and Refs. [23}{26]for alternate perspectives). 

Apart from primordial black holes, the other possible scenario in which small 
black holes might occur is in particle collisions if there are large extra dimensions. 
Large extra dimension scenarios were introduced initially to provide an alternate, 
geometric, resolution of the hierarchy problem. The idea that our four-dimensional 
(4D) Planck scale is derived from a higher dimensional Planck mass close to the 
Standard Model scald248 we then live on a 4D brane embedded in a higher dimen- 
sional spacetime. In such scenarios, it is easier to form black holes in particle col- 
lisions (see e.g. Refs. BIH33). Our relatively high Planck scale, M? = 1/87Gy, is 
then the result of a geometric hierarchy coming from an integration over the extra 
dimensions. 

In a previous paper 24] we computed the probability for seeded decay with a 
brane black hole, following the notion that small black holes could also occur in 
particle collisions if there are large extra dimensions. As in the straightforward 4D 
case, we found the decay rate to be significantly enhanced over the Hawking evapo- 
ration rate for a range of small mass black holes. However, we did not compare the 
seeded nucleation process to an unseeded brane Coleman-de Luccia (CDL) equiva- 
lent rate, hence we could not clarify the extent to which enhancement of tunnelling 
was due to the black hole, or the extra dimensions. In this work, we address this 
question, and explore the vacuum decay of a brane scalar in the absence of any 
black hole. The instantons we will consider will be the true brane equivalents of the 
CDL solution. Early work on brane instanton¢=24 focussed largely on construct- 
ing the Euclidean solutions and presented results on the action within the thin 
wall approximation, subsequent work either focussed on compact instantons in, or 
near, the thin wall limit BI or approximate Hawking—Moss-type instantong28l and 
bulk scalar instantons B? See also some interesting ideas on Randall-Sundrum (RS) 
brane decay by 5D “bubble of nothing” type processes 4044) as well as instantons in 
Dvali-Gabadadze-Porrati (DGP) #2) 

In this paper, we consider vacuum decay of a scalar field localised on a brane 
embedded in a 5D anti-de Sitter (AdS) bulk. We first review the derivation of the 
instanton equations in Sec. P| then present numerical solutions for the scalar CDL- 
equivalent brane instanton in Sec. B] In Sec. [4] we turn to a computation of the 
action, showing how to renormalise the instanton action properly, and computing 
the action for a range of potentials and Planck mass hierarchies before concluding 
in Sec. D] 
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2. The Instanton Equations of Motion 


In the RS model, spacetime is 5D with a negative cosmological constant living in 
the bulk. This negative curvature of spacetime causes a localisation of the graviton 
on the brane, the background solution being a brane with energy and tension equal 
and precisely tuned to the bulk cosmological constant, giving a flat brane at z = 0: 


ds? = e?n da" dx” — dz”, (2.1) 


where £? = —6/As; is the AdS curvature scale and 7,,, the Minkowski metric. The 
local negative curvature of the bulk supports the brane tension ø that is easily 
calculated from the Israel junction condition 


1 3 
KY = gw = Kis — KO Nw = pee = ATG5 Nuw. (2.2) 


One can add energy momentum to the brane, for example a “brane cosmological 
constant,” so that ø is greater or less than the critical value BI a, cosmological 
fluid, or a perturbative localised source. In all cases, the intuitive visualisation of 
brane matter is that it causes the braneworld to bend as first pointed out by Garriga 
and Tanaka”? (see also Refs. [53] and [54). 

We are interested here in pure false vacuum decay, i.e. the brane equivalent 
of a CDL instanton that is a Euclidean solution to the Einstein plus brane scalar 
field equations that has O(4) symmetry on the brane. This level of symmetry is 
mathematically equivalent to a cosmological braneworld solution: there is a brane 
coordinate T upon which the brane solution depends, and a coordinate that tracks 
the warping in the bulk. If we assume that the full brane plus bulk solution also 
has O(4) symmetry, then a “generalised Birkhoff theorem” applies™ and the bulk 
equations of motion can be fully integrated with the brane now following a tra- 
jectory in the bulk consistent with the local energy-momentum of the instanton 
solution (for proof see Refs. [35] [36] and [50). 

To find these equations of motion, we take a simple scalar field Lagrangian on 
the brane in the Wick-rotated Euclidean signature geometry: 


Lg = Thupu +V(¢). (2.3) 


The general bulk admitting an O(4) symmetric brane solution is a Schwarzschild- 
AdS black hole B550 however, as we are computing the brane equivalent of the CDL 
instanton, we will take the pure AdSs spacetime in the bulk 


d 2 2 
T Eed he)=1+3 


sian = h(r)dt? + WA B 


(2.4) 


as a bulk black hole induces a cosmological radiation source on the brane 2951) 
The brane traces out a submanifold in (2.4) that can be parametrised by intrinsic 
coordinates {7, 0°} (a = 1,2,3): 


X” = (t(r), a(T), 0%), (2.5) 
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where 7 is chosen to be the proper time parameter on the brane and a(r) the radial 
trajectory of the brane. 


"E 
h? + —=1, (2.6) 


so that the induced brane metric is identical to the CDL geometry: 
ds? ane = aT + a?(7)dN?,,. (2.7) 


The scalar field depends only on 7, and the energy-momentum is readily found to 
be 


l. 3E 
Tyr = Ve] = Ta 
TFV ae = Te 
i aa (2.8) 
Tag = = h2 a = ——— Jaé 
6 lo+vs 5? lo B = FG, 908 
that sources the brane trajectory. 
The Israel junction equations are then 
1 h! 
Kt = t (a- 2M) =% ~ 27, 
ht 2 
(2.9) 


th 
Kye = 7 Ga8 = —EGap, 


usually expressed in the cosmological format of Friedmann and conservation of 
energy-momentum equations: 


alle, 
a a L 
(2.10) 
E He- T). 


For numerical integration of the scalar P it is more useful to use the Ray- 
chaudhuri equation, and substituting in the form of the energy-momentum (2.8) 
we finally arrive at the full set of brane scalar instanton equations: 


a\> 1 8nGy des AnrG yl \? ba 
Te a 


-= SGN y 4 pp (z (v- za) (v +80"), (2.11) 


3 3 
~ 3a; OV 
ot T? = 36° 

where we have substituted the Newton constant Gy = G;/¢ in the gravitational 


coupling. These are precisely the Shiromizu-Maeda-Sasaki#= equations with van- 
ishing Weyl term, also analysed in Ref. [38] for the Hawking—Moss case. As £ drops, 
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gravity becomes more strongly localised on the brane, hence the 4D limit is £ — 0, 
and (2.11) become the 4D instanton equations. 

It is also worth noting that the critical RS brane (with V = ¢ = 0) has å = 1. 
This leads to the brane trajectory 


r=a(r)=7, t(r)= “log (: 4 =) (2.12) 


in terms of the original coordinates (2.4). This is a less familiar form for the critical 
RS brane, obtained because we are solving for the brane in bulk global coordinates, 
rather than the usual Poincaré patch. The trajectory can easily be transformed to 
its familiar form using 


e” = —— at = el rnh, (2.13) 
ip 


(2 


where n4 is the unit vector in four dimensions. 


3. The Scalar Brane Instanton 


In order to investigate vacuum decay, we use two basic model scalar potentials. 
The first is a standard quartic potential V}, with a potential barrier between a false 
and true vacuum. It is convenient to parametrise this potential with the value of 
ġ = oy at the maximum and ¢ = ġy at the global minimum: 
H ¢ g 
— — z (¢v + om) + ovom]. (3.1) 
4 3 2 

The potential vanishes at the false vacuum @ = 0 and the value at the true 
vacuum is 


Va) = 9 


Vilov) = Z 4% 2¢m — v). (3.2) 


Note that since we require Vj(¢v) < 0, dv > 2¢ar. 

The second potential we wish to investigate, Vp, more closely approximates the 
Higgs potential. The form of this potential has one local minimum and a barrier, 
where on the far side the potential does not turn up again until it reaches very high 
field values. This allows for the possibility of a phase transition and the nucleation 
of a true vacuum bubble. The potential takes the form 


Va(¢) = = est (d)o", (3.3) 


where the effective coupling 


ks -sf (1 An (1m ir) } (3.4) 


g ~ 1075 is a constant that can be used to tune to the potential to closely fit the 
Standard Model Higgs potential. 
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In each case, we integrate from the centre of the instanton, T = 0, looking 
for a solution that asymptotes the flat critical RS trajectory (2.12). However, note 
that because we set boundary conditions at T = 0 of a= 0, a= 1 and ġ = 0, the 
flat geometry at large T is o — dpy, the value of ¢ in the false vacuum, a — T+c — 
integrating through the bubble wall produces an offset in the value of r relative to t. 
While this is not particularly relevant to the form of the bubble solution, for which 
a(T) is important, it is a crucial observation for the computation of the action, as 
we will return to in Sec. [4] 

The quadratic potential (8.1) is particularly useful for exploring the variation 
from thin to thick bubble walls, and for varying backreaction strengths. To illustrate 
this, we present results for two representative potentials, one giving a strongly 
backreacting thin wall, with parameter values g = 1,@y = My, ọm = 0.4Mp, and 
the other a weakly backreacting thick wall with parameter values g = 1/2,dy = 
M,,om = 0.1Mp; in both cases the Planck scales are M5 = 0.4, Mp = 1, hence 
the bulk AdS lengthscale is 2 = 1/M = 125/8. Figure [I] shows the potential V} 
for these two choices of parameters; note the thin wall potential (shown in blue) 
has a significant potential barrier between the vacua, but less well represents a 
Higgs-type potential, whereas the thick wall potential (shown in red) more closely 
resembles the Higgs potential, having a very small barrier relative to the global 
minimum. 

The scalar field solution is shown in Fig. [2] and demonstrates clearly the distinc- 
tion between the potentials: the thin wall has a clear, sharp transition from false 
to true vacuum around T ~ 25, whereas the thick wall does not even reach the true 
vacuum by the centre of the bubble. The effect of the bubble on the embedding 
of the brane is shown in Fig. [3| The strongly backreacting thin wall brane shows 
the transition between the flat RS critical asymptotic false vacuum brane, and the 
sub-critical true vacuum AdS embedding in the interior of the brane. The weakly 
interacting thick wall has a much less significant displacement, and does not reach 
the spherical shape of the sub-critical brane. 


Val) Val) 
0.0g00s 
0.02 -0.05 0.05 0.10 0.16 020 
0.02 -0.00005 
0.01 -0.00010 
0.01 
= $ 
-0.01 
0.5 1.0 $ 
-0.02 
-0.01 
-0.03 


Fig. 1. (Colour online) The V, potentials referred to in the text. On the left in blue with ọm = 
0.4, and ¢y = 1 (with Mp = 1), corresponding to a well-defined bubble wall. On the right in red 
the potential more closely approximated the Higgs potential, with ¢j, = 0.1, and corresponds to 
a thick wall bubble. 
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10 20 30 40 50 10 
Fig. 2. 


20 30 40 50 g 


(Colour online) The scalar field solution for the potentials shown in Fig. [] Once again, 


blue corresponds to the thin wall bubble, here clearly seen as a step in ¢, and red to the thick 


wall bubble. 
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Fig. 3. 


z 
15 


(Colour online) The geometry of the brane with bubble embedding shown in Poincaré 


coordinates, as is usual for the flat RS brane. (Colour scheme as Fig. [I]) 


4. Computation of the Action 


Having found the Euclidean brane bubble solutions, we now need to compute their 
action, in order to find the leading order exponential behaviour of the tunnelling 
probability. We first derive the action for a general brane scalar solution, given 
a large r cutoff, then discuss the proper background subtraction. The Euclidean 


action is given by 
1 2K 
d°a(—Rs + 2A +f d'a] 
8TGs5 s ae 5) amt 87Gs 


dz &? 1 
= —— dz |=- -(V 
fa G5? + fan z 6 3, m 
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that is formally infinite for the background false vacuum critical brane solution. 
Note that this action is written in the Gibbons—Hawking boundary format, with 
the brane being a boundary of a 5D manifold, the normal n, pointing in to the 
manifold — this expression manifestly includes the Zə symmetry of the brane. 

In order to find the instanton action, we first apply a cutoff well outside the 
radius of the bubble. We define the cutoff by a(tz) = R, and bound the bulk 
coordinates by r < R, and t < t,(7R), where ty is the value of t on the brane found by 
integrating the relation (43). Note that t is also bounded below by t,(0). We obtain 


2 tè(R) R TR f2 1 
Sp a at f drr? + an? f dra?(r) a —=(V+o) 
0 0 


= . (4.2) 
Gsl? Ji (0) 3 


Now, whereas the bulk integral is naturally expressed in terms of the bulk coordi- 
nates t and r, the brane integral and the instanton solution are naturally expressed 
in terms of the intrinsic coordinate r. While we can easily identify r = a(r), the 
relation to the bulk time coordinate is differential: 


dt Ea(t 
1+ A 


Using this relation, we can rearrange the bulk integral, integrating first with respect 
to r, then translating to a T integral to finally obtain 


12 TR a? . 
Sr= =| dr ——,[¢? — 2V — 20], (4.4) 
3 0 a 


where TR is defined as a(TR) = R. This integral now is in a simple “brane” format, 
and we can easily insert in the solutions of the scalar instanton equations. The inte- 
gral diverges as 7}, for large R ~ TR, however, outside the bubble, both instanton 
and false vacuum branes are identical, thus once we subtract the background false 
vacuum action this divergence will be removed. 

To subtract the background false vacuum a crucial observation is that the false 
vacuum action is not obtained simply by deleting all but the ø term in (44), since 
not only is a(r) different, but also the value of 7 at which the brane radius becomes 
equal to R (see Fig. Æ). We must therefore perform one final manipulation to get 
the instanton action. The critical false vacuum brane action is 


—2n2 (TR a3(r'\odr’ —2r? [P aoda 
Srv = ip a = | os (4.5) 
3 0 atr) 3 0 a 
1+ 55 1+5 


but now that this is expressed as an integral over a, we can compare this to the a 
integral for the bubble: 


2 pR 3o. 
Sbub = =| L Z $ — 2V — 20] (4.6) 
3 0 a a 


zl 
DEE 
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Fig. 4. (Colour online) The offset between a and r for the thin and thick wall. (Colour scheme 
as Fig. []) 


and using this expression gives the final subtracted action for the instanton as 


Qn? [Ed 3 19 
B=Sp-Spy=— | Z |5-vi(a—vo 
3 Jo à a 2 
I+ a 
2r? TR a? 2 . 

0 
Li 


now expressed as an integral over the brane time-coordinate (and numerical inte- 
gration parameter) r. This action can now be exponentiated to give the dominant 
contribution to the probability of vacuum decay. 

Figure] shows the tunnelling exponent for the potential V} with the parameter 
sets considered in Sec. [3] these are plotted as a function of the mass parameter 
M; = MŽ /3¢-1/3 which determines the strength of gravity in five dimensions. The 


50000 90 
— g=1.0 — 9=0.5 
--- g=2.0 —-- g=0.75 
Peete g=3.0 _| S442 g=1.0 
10000 80 
5000 4 
70 4 
1000 4 
m m 604 
500 5 
504 
100 4 paeer tase SRSeSSeSeSSse= 
aa 40 4 
10 T T T T 30 T T T T 
0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 
Ms/Mp Ms/Mp 


Fig. 5. The vacuum decay exponent B for the quadratic potential plotted as a function of M5 
for barriers with day = 0.4Mp (left) and ġm = 0.1Mp (right). The exponent approaches the 4D 
value as Ms approaches the 4D Planck mass My. 
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0 T T T T 
0.0 0.2 0.4 0.6 0.8 1.0 


Ms/Mp 


Fig. 6. The vacuum decay exponent B plotted as a function of Ms for Higgs potentials with a 
range of values or the instability scale A. There is no dependence on the extra dimension. 


barrier is at ọm = 0.4M, and ġm = 0.1Mp. These test case examples show a 
reduction in B, hence an increase in the vacuum decay rate, due to the increasing 
influence of the extra dimension. 

The edge of the plots denotes a minimum value of Ms beyond which the numer- 
ical solutions cease to exist. Close to this limit, the total surface tension on the 
brane becomes negative near the centre of the bubble. Note that the allowed range 
of M; is narrow, as in the examples plotted above, therefore does not correspond to 
a significant hierarchy. Therefore adding an extra dimension only affects the decay 
rate in very specialised situations. 

We also show the tunnelling exponent for the Higgs-style potential Vaz, with 
parameters chosen within the Standard Model range in Fig. [6] The Higgs potential 
is small at the Planck scale because the parameter g in the potential is small. Conse- 
quently, vacuum decay rates with the Higgs potential show no obvious dependence 
on the extra dimensions. 


5. Summary 


To sum up: we have found instanton solutions for a brane scalar field representing 
vacuum decay from a critical RS flat brane. We explored general bubble solutions, 
as well as an approximate Higgs potential. We calculated the tunnelling exponent 
for a range of warping in the extra dimension, and compared it to that of a phase 
transition in 4D asymptotically flat space. The influence of the fifth dimension on 
tunnelling rates is relatively minor, except for a strongly backreacting bubble. 

A Higgs-style potential was also considered, however, for realistic parameter 
ranges, the impact of the extra dimension was negligible. This is to be contrasted 
to the case of vacuum decay seeded by primordial brane black holes, as in Ref. B4 
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where the decay rate is significant. We conclude that, rather like the 4D case, black 
holes are required to produce significant decay rates. 

One interesting feature of our numerical solutions was that they had a sharp 
cut-off in the allowed value of M5, due to the brane tension becoming negative. This 
is possibly due to the fact we integrate out from T = a = 0, hence this method does 
not allow for a wormhole-type solution where the brane transitions from positive 
to negative tension as in Refs. [40] and [41] It might be interesting to consider this 
further. 
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